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We investigate the formation of caustics in the Dirac-Born-Infeld type scalar field systems for generic classes 
of potentials, viz., massive rolling scalar with potential, V{4>) = Vbe^a and inverse power-law potentials 
with V{<l)) — Vo/4i" , < n < 2. It is found that in the case of exponentially decreasing rolling massive scalar 
field potential, there are multi-valued regions and regions of likely to be caustics in the field configuration. 
However there are no caustics in the case of exponentially increasing potential. We show that the formation of 
caustics is inevitable for the inverse power-law potentials under consideration in Minkowski space time whereas 
caustics do not form in this case in the FRW universe. 

I. INTRODUCTION 

The discovery of the late time accelerated expansion of universe is one of the most surprising findings in modern cosmologyiH 
and posses a serious challenge to fundamental theories of physics. The resolution of cosmic acceleration riddle similar to the 
phenomenon of black body radiation might unveil new secrets of nature. 

It is now widely believed that the late time acceleration of universe is due to an exotic form of energy with large negative 
pressure known as dark energy which is the dominant fraction of the energy content of present universe ||3-Sl- The simplest 
dark energy model based upon cosmological constant is faced with fine tuning problem of an acceptable level. As an alternative 
to cosmological constant, a variety of scalar field models is proposed in recent years to provide a viable explanation for the phe- 
nomenon of late time cosmic acceleration IH-l^. Though the scalar field models have limited predictive power but nevertheless 
can be of interest in case they can exhibit some generic features allowing to alleviate the fine tuning and coincidence problems or 
could be motivated from a fundamental theory of high energy physics. The Dirac-Born-Infeld (DBI) scalar field model is string 
inspired and certainly invites attention ifioll . Unlike generic quintessence models with tracker solutions, there exists no solution 
which can mimic scaling matter/radiation regime in case of the tachyon field ifTlUigll . These models necessarily belong to the 
class of thawing models: At early times, the expansion dynamics is governed by the background fluid whereas the tachyon field 
remains subdominant and frozen. However, as the background energy density redshifts and becomes comparable to the field 
energy density, the field begins to evolve and subsequently overtakes the background to become the dominant component of 
universe. 

Tachyon models do admit scaling solution in presence of a hypothetical barotropic fluid with negative equation of state. 
Tachyon fields can be classified by the asymptotic behavior of their potentials for large values of the field; (i) 1^(0) faster 
than 1/0^ for — > od. In this case dark matter like solution is a late time attractor. Dark energy may arise in this case as 
a transient phenomenon, (ii) V{(p) — ?> slower than l/cfy^ for — > oo ; these models give rise to dark energy as late time 
attractor. The two classes are separated by V{(l>) ^ 1/4P' which is a scaling potential with w{4)) = constant. The present state 
of observations does not allow to distinguish amongst various scalar field models and leaves the dark energy metamorphosis as 
a future challenge for observational cosmology. 

The formation of caustics in field profile in the mass free space is an undesirable consequence in the field theoretical models 
in cosmology as it indicates the failure of physical theories to explain the evolution of field in that particular region. Thus the 
study of formation of caustics in the field configuration is one of the best techniques to investigate the fundamental shortcoming 
of the field theory for a specific potential. 

Inspite of the exiting features of cosmological dynamics based upon DBI scalars, it might happen that these models lead 
to formation of caustics where the second and the higher-order derivatives of the field become singular As demonstrated in 
Ref. lITsll . caustics inevitably form in tachyon system with potentials decaying faster than (f)^"^ at infinity. We do not know 
whether caustics are generic prediction of string theory or appear as a result of the derivative truncation leading to the DBI 
action. It remains to extend the analysis of Ref. [15,1 to the case of inverse power law potentials, V ^ 0^" with < n < 2 
analysed in Ref. ll22ll . Caustics normally form in systems with pressureless dust which is mimicked by tachyon field with run 
away potentials. It is therefore quite likely that caustics may not develop in Born-Infeld systems with a ground state at a finite 
value of the field. The rolling massive scalar potential, V{4)) = Vbe 2 belongs to this category. 

In this paper, we address the issue of caustic formation for DBI scalar field systems for inverse power potentials giving rise 
to dark energy as a late time attractor and the massive rolling scalar In the section II, we present the DBI type general scalar 
field equation for the \ + \ flat space time and review the formalism of field dynamics related to caustic formation based upon 
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Ref. ifTsll . We also extend the analysis to the case of isotopic and homogeneous expanding universe to see the formation of 
caustics in the field in more real situation. In section III, we discuss our numerical results. To have a comprehensive view 
of the scalar field pattern, the general field equations are solved numerically as well analytically in the homogeneous and 
inhomogeneous situations for the two different classes of field potentials mentioned above. For the inhomogeneous case, we 
solve the evolution equations analytically using the method of characteristics |fl5| to check the formation of caustics in the field 
profile using the geometrical interpretation of the solution. Main results of our analysis are summarized in the last section. 




FIG. 1 : The numerical solution of the equation ( I26t (red line) and the analytic solution of equation ( I27t (black line) for different values of 
A'l and (j>o- The solid line is for the initial value of field (j>o = 0.001 and the dotted line is for (j>o = 0.1. Three different sets are due to three 
different values of AI = 1, 0.5 and 0.1 along the increasing values of t. The discrepancy between these two solutions is obvious after some 
initial period depending upon the values of AI and (po- 



ll. THE DIRAC-BORN-INFELD TYPE SCALAR FIELD SYSTEM 
The DBI type action for a scalar field 0, referred to as tachyon field hereafter, is given by iH [T 



5 = -y d4^y^y(0)v/TT9;:^, cd 

where V{(j>) is the potential of the field cj). The field equation derived from the action ([T]) is, 

where the covariant derivative of the field (j> with respect to the metric g^j/ is denoted by V^. 

A. Minkowski 1+1 dimensional analysis 

In what follows, we shall first consider the field equations in 1 + 1 dimensional Minkowski space. Caustics formation can 
easily be analysed in this case. If caustics do not form in Minkowski space time, we would expect the result to hold in the 
expanding universe as expansion should work against the formation of caustics. However, in the opposite case, it essential to 
incorporate expansion to reach the final conclusion. 

In this case, the equation ^ takes the form. 



= 1 



(3) 



Here the time and space derivatives of the field </> are indicated by the dot and prime over </> respectively. 

In order to understand the time and space evolution of the field (/> in 1 + 1 dimension under a specific field potential, let us first 
consider some general aspects related to the equation (O for the homogeneous and inhomogeneous fields. 

If the scalar field is homogeneous, then 0' = 0, and hence the equation (O for such a field simplifies to, 

4> VM) 



1 - 02 V{<j)) 
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For a particular potential, the time evolution of the homogeneous field can be obtained from this equation In the case of 
an inhomogeneous scalar field, <f>', (j>" > 0, however they should be sufficiently less then unity for a realistic scalar field in 
cosmology [15]. Under this condition, the equation ([3]) can be written as. 



= 1- 



1/(0) 



^(0)Q(</'), 



(5) 



where. 



P{<j>) = 1 - <^ 

Q{4>) = 0" - 



V{ct>) ■ 



(6) 



The equation of P{(t>) contains both the time and space derivatives of the field and therefore this parameter can be used as 
an indicator of the pattern of evolution of the field governed by the particular field potential. For example, if the field 4>{x, t) 
rapidly approaches the configuration P((/)) = 0, it then indicates that the time and space evolution of the field is such that it 
rapidly approaches 0^ — 0'^ = 1. Alternatively, we may say that, under this situation the field (j>{x, t) is almost similar to some 
subsidiary field 9{x, t) whose time and space evolution is constrained by the equation ifTsll . 



1. 



(7) 



The field may approaches to the configuration ?{(()) = from both directions, viz. from above or below zero, depending upon 
the initial state of the field. If the initial state of the field is such that initially P(0) > 0, then P{(t>) will asymptotically approach 
to zero from above, otherwise it will tend to be zero from below lITsll . Similarly, if P{(j)) ~ 1, it implies that, 4? w (f)'"^, i.e., the 
type of evolution of the field 0(.t, t) with respect to space and time are nearly equal. In this case also, we may consider that the 
(/)(a;, t) is almost similar to some subsidiary field t) whose time and space variations are related by the equation. 



■ 2 12 



= 0. 



(8) 



The extensive numerical solutions of the equation (|5]l for the scalar field with the two different classes of field potentials of our 
interest clearly showed the above behavior of P{(t>)- Thus the equation Q has two attractors, one for P{(t>) — and the other 
for P{4)) — 1, depending on the field potential and hence P{(t>) can be used to define the following relation, 



= ^l + 0'2_P(0). (9) 
If we consider that the scalar field 4'{x,t) itself satisfies the subsidiary field equations (|7]l and ^ under different conditions. 




FIG. 2: The numerical solutions of the equation i29\ for different values of M and 4>o- Black, red and green lines are for M = 0.1, 0.5 and 1 
respectively. The solid line represents (j>Q =0 and the dotted line = 0.1. For all these plots we consider 0' = 0.01 and (j>" = 0.02. The blue 
line represents the solution with (/>o = 0, M = 0.5, 0' = 0.002 and = 0.001. For all cases time variation of the field is very slow initially and 
then it increases rapidly as time passages. 



then analytic solutions of these two equations correspond to free relativistic massive wave propagation along the characteristics 
of the field (j){x, t). From the particle point of view, equations Q and (|8]l can be viewed as the dynamical descriptions of the 
motion of free massive particle under different situations. From this point of view, the trajectories x{t, q), (where q is the initial 
spatial coordinate of the field, i. e. x(0, q) ~ q) of individual particles and the corresponding evolution of the field configuration 



FIG. 3: The variation patterns of P(^) and Q{(t>) witli respect to time for tlie potential i25l corresponding to tfie initial conditions of the figure 
|2] The falling lines are for P{(j}) and raising lines are for Q{(t>)- The overall pattern are same for all initial conditions. 



can be obtain from the parameterized solutions of the characteristic equations of two variables, viz., q and the affine parameter 
s along each curve, which are given by 115,1 , 

dx 

— (g,s)=P,0' = 20'; xiq,0) = q. (10) 

J(g,s) = P,0 = -20; i(g,0) = 0. (11) 
as 

^(g,s) = 0'P,(/.' + 0P,(/. = 20'2-202. ^((^,o) = 0,(9). (12) 

Solving these equations and using equation (|9]l we may write, 

x{q,s) = q + 2s(f>i^q. (13) 



t(g,s) = -2,s^l + 02^-P(0). (14) 

<Piq, s) = Ml), for P(0) = 1 and 0(<z, .s) - ^,{q) - 2s, for P(</)) - 0. (15) 

Where (j)i^q = <p'{q, s). Eliminating the affine parameter s from the above equations, we obtain the parametric solution for the 
field (j){x, t) as, 

xiq,t)=q- '^^''^ t. (16) 

<j){q,t) = (P,{q), for P{(f>) = 1 and (f>{q, t) = </),(g) + , for P(</.) 0. (17) 

The equation ( fT6l ) will provide us the trajectories of the free massive particles in the field and the equation ( fTTb will give us the 
pattern of evolution of the field with respect to time fisll . Hence these two equations will give us the geometrical form to see 
the formation of caustics in the field profile and the evolution of the field configuration with any presumed initial field profile 
under a specific field potential. In what follows, we shall consider the potentials of our interest for the specific solutions of the 
equations (|4|i and (|5]|. Although our main interest is the solution of equation for the completeness of the picture we present 
the solution of equation (|4|i also. From the solutions of equation (|5]) for the following three different field potentials, we shall 
obtain the characteristics curves from equation (fT6b which will show clearly whether there are caustics in the field configuration. 
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B. Generalization to the case of expanding universe 

So far, for simplicity we have considered the field in 1 + 1 dimensional Minkowski space time. However for the real cos- 
mological applications, it is necessary to consider the field in 3 + 1 dimensions in the expanding universe. It should be noted 
that the extension of the field from 1 + 1 dimension to the 3+1 dimensions is a matter of analogical enhancement of the one 
dimensional spatial coordinate to the three dimensional spatial coordinates. Hence the basic form of all equations discussed 
above will remain same in 3 + 1 dimensions, also except that the space derivative has to be replaced by the nabla operator V, 
and the spatial coordinates x and q by the vectors r and q respectively in the relevant equations. Similarly if we take a spatially 
flat Friedmann-Robertson- Walker (FRW) metric with a scale factor a(t), then equation for the scalar field in an expanding 
universe can be obtain from the equation (|2|i as, 

'<h ■ VM) 

3i/0+-p^ = O, (18) 



1-02 ^ 

where H = a/ais the Hubble rate and can be expressed in the scalar field dominated expanding universe as. 



o SnG V((h) 

= — (19) 



The equation ( flSl l is equivalent to the equation (|4]i of 1 + 1 dimension. If some small inhomogeneous perturbations arise in the 
scalar field of the expanding universe, then these perturbations will generate the small perturbations in the FRW metric which 
can be expressed by the Newtonian gravitational potential $0(^7 r) as |[l5il23ll . 

ds^ = -(1 + 2^G)dt^ + (1 - 2<^G)a^{t)dr^, (20) 

where we have use the metric convention as (— , +, +, +). If for a particular potential of the scalar field of expanding universe 
with small inhomogeneous field perturbations lead to the state of the field as in the cases of the equations (jTji and (O of the 1 + 1 
dimension, then the corresponding equations for the field itself of the expanding universe represented by the metric equation 
( I20I 1 are given by IH], 

</.2-l(V0)2 = H-2$G- (21) 

02 - \{V<l>f = 0. (22) 

If we solve these two equations by the method of characteristics as in the case of 1 + 1 dimension the corresponding equations 
of (HSll and dnli would be, 

r(q,i) =^.--k—p '^^q'^)^ forP(0) = andr(q,t) = q - -, for P(0) 1. (23) 



l + ^(Vq0)2+2a>G 

(q, t) = 0,(q) + , (l + 2^G)t p^^^ ^ Q ^^^^ ^ ^^^^^ p^^^ ^ ^24) 



The scale factor a{t) assumes different forms at different stages of the evolution of universe. For example, in radiation dominated 
stage, a{t) = aot^^^ and in the scalar field dominated stage a{t) = aot^/^ lfl5il23ll if the tachyon potential vanishes at infinity 
faster than In case of rolling massive scalar or the inverse power-law potentials V{(j)) ^ 1/0" with < n < 2 , dark 

energy is a late time attractor and the scale factor takes the form a{t) = aot2/(3(i+«')) with suitable negative values of w. The 
fluctuations in the tachyon field grow lineally with time while the metric fluctuations remain constant for potentials decreasing 
faster than 1/0^ at infinitv lll5n23il . However, in case of rolling massive scalar and inverse power- law potentials with < n < 2 
that would be of interest to us, both the fluctuations do not grow (jgl IIU 



III. CAUSTIC FORMATION 



In this section we apply the above formalism to investigate the possibility of caustic formation in the DBI tachyon system 
with two generic classes of potentials and present the results of numerical simulation. 
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A. Exponentially decreasing rolling massive potential 
The exponentially decreasing rolling massive scalar field potential is given bv ll22ll . 

F(0) = Foe-^*^'^', (25) 
where Vq and M are constants. For this potential the homogeneous scalar field equation (01 becomes, 

-X- - = 0. (26) 

The numerical solution of this equation is shown in the figure[T|for different values of the arbitrary constant M and initial values 
of field 00. The values of M chosen in this figure are 0.1, 0.5 and 1. On the other hand two values of 0o are considered for this 
figure, which are 0.001 and 0.1. If the time variation of the scalar field is very slow, then we may consider that << 1. In this 
case, equation ( |26] | further simplifies to, 

- 1/20 = 0. (27) 
The solution of this equation is trivial and can be written as, 

= cj)ocosh{Mt), (28) 

where 0o is the constant initial field value. The plots of this equation are also shown in the figure[T]for same values of M and 
00 as for the above case to compare with the numerical solution. The red line indicates the numerical solution and black line 
indicates the analytic one. The solid line is for the initial field 0o = 0.001 and the dotted fine is for 0o = 0.1. Three different sets 
are due to three different values of M = 1, 0.5 and 0.1 along the increasing values of t. It is observed from the figure [T] that the 
assumption of very slow variation of the field with respect to time is correct only for some initial period that depends upon the 
values of M and 0o. For low values M and 0o this period is longer than their higher values. As the time increases beyond this 
period, the discrepancy between the solutions of equations ( |26] | and ( l27b increases rapidly. The overall trend of time variation 
of the field with different possible initial conditions is almost same, the only difference is that, this variation is slower for some 
initial period if the field is started with smaller values of M and 0o than their higher values. 
The inhomogeneous scalar field equation© for this potential can be written as, 

- (l - 0^ + 0'^) (0" + Af 20) . (29) 

The numerical solutions of the equation ( |29] ) for the same values of AI as in the previous cases, but for 0o = and 0. 1 are shown 
in the figure |2l In the figure the black line is for M = 0.1, red line is for M = 0.5 and green line for i\/ = 1. The solid line 
represents for 0o = and the dotted line for 0o = 0. 1 . Since the solution of the inhomogeneous equation ( l29b depends upon the 
values of 0' and 0", so for these plots in this figure we consider 0' = 0.01 and 0" = 0.02. Moreover to check the dependency 
of the solution on 0' and 0", we consider another set of values of these two parameters as 0.002 and 0.001 respectively together 
with M = 0.5 and 0o = 0.1. The result of this solution is shown in this figure by the blue line. We observe that for all cases the 
variation or the evolution of the field with time is very slow in the initial period based on the values of the Af and 0o, then it 
develops rapidly as time passages. For low values of M and 0o this initial period enhances slightly and the time variation is also 
relatively slower for such values. The value of A/ has significant effect over the variation of the field with time than 0o which 
is also noticed for the homogeneous case. Since this is an arbitrary parameter, we can assume any value of it for a solution, 
however its actual value can not be substantially small because for such a small value the time evolution of the field will halt for 
a considerable period of time. The 0' and 0" does not have any effect on the time evolution of the field, whereas they only effect 
to shift the magnitude of field in the same direction of their magnitude variations (i.e. along the increasing or decreasing order). 
Thus for our remaining discussion, related to this field potential, we will not consider again the effect 0' and 0" over the field. 

The variations of P(0) and (3(0) with respect to time that are obtained from the numerical solutions for 4){x, t) are plotted in 
the figure[3] In this figure, the falling lines are for P(0) and raising are for (5(0). The color and style of the line represent the 
corresponding values of M and 0o as in the case of the figure |2l As mentioned earlier the state of the field with respect to time 
and space is indicated by the behavior of the P(0) at the corresponding time, so, here we are mainly interested in P(0) only. 

From the pattern of variation of P(0) with time we observed that during some initial period, which depends upon the value 
of A/ and 0o, the value of P(0) is w 1 and beyond this period its value falls rapidly to zero. The length of this initial period is 
more for lower values of M and 0o (effect is prominent for the parameter M as mentioned earlier), otherwise the behaviour of 
time variation of P(0) is similar for all initial conditions. To be more specific, we consider the solution of P(0) with A/ = 1 and 
00 = 0, and the initial field configuration as 0i('?) = exp{—q^) ifTsll . then we obtain the characteristic curves from equation ( fTSl ) 
which are shown in the left panel of the figured From these characteristic curves it is interesting to note that there are likely to 
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FIG. 4: The particle trajectories in 1 + 1 Minkowsld space obtain from equation J16l l (left) and in FRW expanding universe obtain from 
equation l |23t (right) with potential l |25t and initial field profile 0i(q) = exp{—q^). Both panels of figure are obtain with the initial condition 
M = 1 and (t>Q = 0. 



be caustics as well as multi-valued regions in the field profile. Apart from these, there are regions of twisting of the characteristic 
curves. The regions of caustics and multi-valued start at the points x = ±3. As the behaviour of the P(0) is similar for all 
initial conditions mentioned above, therefore we will get the similar characteristic curves with caustics and multi-valued regions 
at different locations corresponding to the values of M and 0o. Obviously for smaller values of M and 0o these location will 
shift to the higher value of t. But it is sure that there is no way to avoid these unphysical regions for any value of initial field 
parameters. Thus we may infer that the formation of caustics and multi-valued regions are independent of the initial conditions 
of the field. 

The evolution of the field configuration (p{x, t) with (f)i{q) = exp{—q^) and for the same initial condition as in the case of 
characteristic curves are shown in the figure |5] at different times. It should be pointed out that to draw the field configuration 
at t = 5 we consider that P{4>) = (which is almost the case. Beyond this value of t, P{(f>) is obviously zero). Within the 
period, when P{4)) falls from almost one to zero, field evolves very rapidly with time, which is not significant to be visualized 
geometrically. The same observation can be made for the evolution of the field configuration with different initial conditions and 
consequently with the different values of time after which P(0) falls to zero. 

Heuristically speaking, expansion works against caustics formation. It is therefore necessary to investigate the corresponding 
situation in case of the expanding universe using equations ( |23] l and (|24] |. Without going into detail of the scalar field dynamics 
in the expanding universe under this potential, we assume that the field behaves in a manner equivalent to 1 H- 1 dimensional case 
(which is specified by the behavior of Since the average matter density of the universe is quite low to obtain a substantial 

value of the Newtonian gravitational potential $g(^i r) [24, :25J of the universe, we shall neglect in comparison to unity. 
Secondly, as dust like solution is a late attractor in the present case, we assume a{t) ^ t^/'^. Under these considerations, we 
obtain the characteristic curves from equation (l23T l with the initial field profile as mentioned above are shown in the right panel 
of the figure |4] for the one dimensional space. It is observed that the patterns of the characteristics curves are different in this 
case from the 1 H- 1 dimensional Minkowski space time and the caustics are more distinctly formed in the field profile if the 
field behave in a similar manner as in the 1 H- 1 dimension with this exponentially decreasing rolling massive potential. From 
the numerical calculations under the above considerations we found that the evolution patterns of the field obtained from the 
equation ( |24] | remain nearly same as its initial profile. We should note that the field rolls slowly for a short while near the origin 
but quickly enters the fast roll to mimic dark matter like regime described by a{t) ~ t^/^. Using equations (18) and (19), we 
have verified that situation depicted in figure |4]does not change if the detailed field dynamics is taken into account. 
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FIG. 5: The evolution of the field configuration (^(i, t) with potential i25l and (l3i{q) = exp{—q^) at different times for the initial condition 
M = 1 and 0o = 0. 



B. Exponentially increasing rolling massive potential 
Now we consider rolling massive scalar potential V{(f>) given by ll20ll22ll . 

= Foe^^^'*^ (30) 
In this case, the homogeneous field equation can be written as, 

^. +M^(P = 0. (31) 

The numerical solutions of this equation is shown in the figure |6] by dotted line with three different colors for three different 
values of M as in case of previous potential, viz., black for 0. 1 , red for 0.5 and green for 1 . The upper set is for 00 = 1 and lower 
set is for (f>o = 0.5. The field <f> for this potential is oscillatory in nature. Since (p is oscillatory just like a simple harmonic motion 
and its obvious equation can be obtain by considering << 1, which is in this case is, 

4> + M^<f) = 0, (32) 

and the nature of is the solution of this equation, 



<j) = (j)ocos{Mt), 



(33) 
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FIG. 6: The numerical solutions of the equation Olt (dotted line) and the solutions of equation l l32b (solid line) for three different values of M 
= 0.1, 0.5 and 1, which are represented by black, red and green colors respectively. The upper set is for (j)o = 1 and lower set is for 0o = 0.5. 
There is a good agreement between the numerical and analytic solutions although there is gradual phase shift between them as time passages 
depending on the values of M and No- 
where 00 is a constant, the field amplitude. The solution ( l33T l are presented in the figure |6] for same values of M and 4>o as 
described above. It is seen that the solutions of equations dSTT i and (|32] | are in good agreement, however there is gradual phase 
shift between them because the second solution is due to slow variation of 0. The rate of this phase shift depends upon the values 
of M and (j>Q, which is slower for the smaller values of these two initial field parameters. On the other hand, the time period 
of oscillation of the field depends on these two parameters in the opposite way, i. e. the time period increases with decreasing 
values of them. As in the case of the previous potential, in this case also the effect of M on the field is more prominent than (pQ. 
The amplitude of oscillation solely depends on the value of (f>Q as it is clear from the figure. 
The inhomogeneous scalar field equation ^ in this case can be written as. 



The numerical solutions of the equation ( |34] | with cj)' ~ 0.01, 0" — 0.02 and with three different values of AI same as above, 
and 00 = 0, 0. 1 are shown in the figure]?] The three colors corresponds to the respective values of M as mentioned above for the 
homogeneous case. Here the solid line indicates the solution for (po = 0, whereas the dotted line for the 0o = 0. 1 . To understand 
the effect of (f>' and <f>" on the inhomogeneous field for this potential we consider the solution of the equation ( [34l ) with the set 
of initial field values as (f>' = 0.002, 0" = 0.001, M = 0.1 and (po = 0, which is shown by the blue line in the figure. We have 
seen that the inhomogeneous field is also oscillatory in nature, time period of which depends upon the values of M and (po in the 
same way as in the homogeneous case, however the small value of M effects significantly on both time period and amplitude 
of oscillation. It should be noted that, the amplitude decreases for higher value of (po as M decreases and there is a phase shift 
between oscillation patterns of field for different values of <j>o, which decreases with decreasing value of M. The effect of (p' 
and 0" on the field is only on the amplitude of its oscillation and we have seen that the amplitude decreases substantially for 
the lower value of space variation of the field. So the role of space variation factor is significant for this potential to avoid the 
unwanted large fluctuation or oscillation of the field for lower value of AI. It also implies that in reality the space variation of 
the field should be sufficiently less that unity as previously mentioned. 

The time variations of P(0) and Q{(f>) corresponding to above solutions of equation ( [34l i are shown in the figure[8] The values 
of P(0) oscillates with time near unity from below with different amplitudes corresponding to the situations described above for 
the solutions of equation ( l34b . For the value of M = 0.1 the amplitude of oscillation of values of P(0) is noticeably away from 
unity in comparison to other higher values of AI if we consider the value of 0' = 0.01 and 0" = 0.02. So for further lower values 
of M the amplitude of oscillation should move far away from unity. On the other hand if we consider the value 0' = 0.002 and 
0" = 0.001 for the same value of M (i.e. 0.1), the oscillation of the values of P(0) is negligible and its values are almost unity 
for all times. But intuitively enough, we may argue that the value of space variation of the field can not be very small beyond 
some limit as that will nullify the inhomogeneity condition of the field and similarly the value AI also can not be very small 
because then field will not be sustainable one for very large amplitude of oscillation. Hence for the real situation the value of M 
and the space variation of the field should be such that oscillation of the field must be negligible. Under such considerations and 
for our present values of M we may take on average the value of P(0) as unity for all times. For the initial field configuration 
= (7^) as in the previous cases and P(0) = 1, the characteristic curves for equation ( fTSI l are shown in the left panel 

of the figure |9l It is clear from the panel of the figure that there are no caustics and there are no multivalued regions in the 
field configuration. This happens due to the almost steady nature of the field (because the average value of the field with time 




(34) 
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FIG. 7: The numerical solutions of the equation ([34} for different values of AI, (j)' and 0". The line colors black, red and green are for M 
= 0.1, 0.5 and 1 respectively. Solid line refer to 00 = and dotted to c^o = 0.1. All these plots are drawn for solutions with (f}' = 0.01 and 0" = 
0.02. The blue line indicates the solution with M = 0.1, 0o =0,(t>' = 0.002 and 0" = 0.001. 




FIG. 8: The time variation patterns of P((^) and Q{<j)) for the potential i30l corresponding to the solutions of equation l l34t shown in the figure 
III 



is a constant as clear from our argument together with the figure |7ll and therefore the evolution pattern of the field for all times 
in this case will be same as in the previous case for i = 0, i. e. same as the initial configuration of the field. As expansion 
works against caustics formation, we should not encounter caustics if expansion of universe is incorporated. Indeed, we have 
numerically checked that there are no caustics in the field profile in the expanding universe in this case, see the right panel of the 
figure |9] 
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FIG. 9: The particle trajectories in 1 + 1 Minkowski space obtain from equation J16b (left) and in FRW expanding universe obtain from 
equation i23\ (right) with potential l |30| l and 4>i{q) = exp{—q^). The figure in the right panel is obtained from the numerical solutions of 
equations dlSt and M9\ of FRW expanding universe with the given potential for one spatial dimension taking constants k — 8ttG, Vq, M as 
unity and the initial field value cj>o = QA. 



C. Inverse power-law potentials 



Finally we consider the inverse power-law potential of the form, given by Il20ll22ll . 

= ^00^", < ??. < 2. (35) 



The late time accelerated expansion of the universe corresponds to < n < 2 112 IL 12211 . and accordingly we restrict the value of 
the exponent 7i within this range for our work. 

We shall first restrict our attention to 1+1 dimensional Minkowski space time. In this case, the homogeneous field equation 
(|4|i is given by, 

^ n 

- - - = 0. (36) 



1-./.2 

The numerical solutions of this equation is shown in the left panel of the figure[TO]for n = 0.5, 1.0 and 1.5 with (po = 1. The 
difference of the solutions decreases for the higher values of n as it is clear from the figure. The inhomogeneous scalar field 
equation (O in this case will take the form, 

?i=(l-0' + 0'')(0" + ^). (37) 

The numerical solutions of the equation (|37] | with cj)' ~ 0.01, (j)" = 0.02 and for n = 0.5, 1.0, 1.5 with 0o = 1 are shown in the 
middle panel of the figure[TO] It should be noted that there is no significant difference between the solutions of homogeneous and 
inhomogeneous field equations for this inverse power-law potential, because the time variation of the field is more prominent 
than the space variation we have considered, which indeed should be for a mass free space. Moreover the field evolution pattern 
is almost independent of the initial condition of the field. However it has little impact on the magnitude of the field as time 
passage in a arbitrary manner as it is seen from the right panel of this figure. So the value of (po would be immaterial for the 
testing of the caustic formation in the field for the inverse power-law potentials. However we take care of the effect of 4>o on 
the variation P(0) to see the way of rolling of its values with time which will be clear from the following discussion. The time 
variations of P(0) and Q{(j>) corresponding to the solutions of 4>{x, t) as shown in the middle and the right panel of the figure 
[To] are shown in the figure [TT] 

From the figure[TT]we observed that values P{(p) for all n rapidly falls to zero from above if the initial value of the field (po 
is not very small. It is interesting to note that, when the value of (pa is very small, P{<p) falls rapidly to zero from below as 
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FIG. 10: The numerical solutions of the homogeneous field equation ( 136b (left) and the inhomogeneous field equation ( I37t (middle) with (l>o = 
1 . The difference between the homogeneous and inhomogeneous solutions is insignificant. The right panel shows the three different solutions 
of the inhomogeneous equation for three different values of (j)o and forn = 1. This panel clearly verifies that the time evolution of the field is 
almost independent of the initial value of the field for the inverse power-law potentials. 



indicated by the third panel of the figure[TT] Notwithstanding, for all cases the value of P{4>) rolls to zero very fast (rate is more 
for lower value of (f>o) as time passage and hence considering the initial field configuration as in the previous cases, the plots of 
the characteristic curves for equation (fTST i are shown in the figure[T2]for n = 0.5, 1.0 and 1.5 with (j>o= I- 

The figure [12] vividly shows the formation of caustics and multi-valued regions in this case. It should be noted that the regions 
of caustics are not steady but slowly moves with time as the exponent n of the potential increases. Due to pattern of variation of 
P{4>), we get the similar pattern of variation of field configuration at different times as in the case of exponentially decreasing 
rolling massive potential as shown in the figure |5] 

As noted earlier, expansion works against caustic formation. The effect might become dramatic in case of inverse power-law 
potential under consideration with dark energy as late attractor. Thus the observed caustic formation in case of potential ( l35b 
may quite be the artifact of Minkowski space time ll26ll . In order to reach the final conclusion in this regard, it is essential to 
incorporate the expansion of universe. The numerical solutions of evolution equations ( flSl ) and ( fT9] l are depicted in the figure 
[13] for all potentials, which are obtained by taking constants k = SttG, Vq, M, n as unity and the initial field value <f>Q = 0.1. 
The left panel provide us the evolution history of the scale factor a{t) and right panel shows the variation of the equation state 
parameter w of the field with time for all three potentials. It is clearly observed that, the time evolution of the scale factor is 
faster in the case of exponentially increasing massive potential than the inverse power-law potential. On the other hand in case 
of exponentially decreasing massive potential the scale factor becomes infinity after some initial period. The field is rolling 
faster to attain steady values in the case of the exponentially increasing and decreasing field potentials. Whereas in the case 
of the inverse power-law potential (0 < n < 2, the plot is shown for n ~ 1) the field is evolving very slowly with time after 
some very brief initial period during which it evolves relatively fast. Though both the models (i. e. exponentially increasing and 
inverse power-law) can account for late time acceleration, the rolling massive scalar requires enormous fine tuning in order to be 
consistent with observation. 

Since dark energy is late time attractor in the present case, we first for simplicity and consistency with the full solution of a{t) 
as mentioned above, assumed that a ^ i'2/3{i+w)^ incorporating all time dark energy effect in the scale factor. Now for the case 
of the scalar field in expanding universe with the inverse power-law potentials under the above mentioned assumptions and with 
this value of the scale factor for w = - 0.9, we obtain the characteristic curves of the field given by the equation (i23] l look like 
as shown in the left panel of the figure [14] It is interesting to note that there are no caustics and multivalued regions in the field 
profile of dark energy dominated universe in contrast to the situation of 1 H- 1 dimensional Minkowski space. Finally at this stage 
we intend to incorporate the dynamics exactly by taking the numerical solutions of equations ( fTSl l and ( [T9] l. The numerical data 
of the scale factor a{t) (see, figure [T3]l can then be used to solve the equation (i23T l to see the formation of caustics in the field 
profile for the given potential. The result of solution of this equation (i23T l, i. e. the characteristic curves obtained by using the 
numerical data of the scale factor is shown in the right panel of the figure [T4]by considering the same initial profile of the field 
as in the previous cases. Both of the panel of this figure is for n = 1. It is apparent from the figure[T4]in the right panel that at 
late times (> 4), the field profile based upon the exact simulation is same as the one (left panel of figure [T4]| that is obtained by 
assuming dark energy domination at all times. This shows that at late times the inverse power-law potentials have exactly same 
effect on the scalar field in the expanding universe as the dark energy has on it. Since the dark energy is a late time effect, so the 
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FIG. 11: The time variation patterns of P(</>) and Q(<j}) corresponding to the solutions shown in the middle and the right panels of the figure 
[Tol 



disagreement of the two panels of this figure during inital evolution is obvious and does not have any significance for us. 

In view of the aforesaid, we consider the models based upon the inverse power-law potentials more viable than the one with 
rolling massive scalar field to explain the present accelerated expansion of our universe. 



IV. CONCLUSION 



In this paper, we have examined the phenomenon of caustic formation in tachyon system with two generic classes of potentials. 
We presented analytical estimates supported by detailed numerical simulation. We show that the time variation of the scalar 
field with exponentially decreasing potential, V — Vbe~ 2 , becomes gradually important as the time elapses for both the 
homogeneous and the inhomogeneous field configurations. There are multi-valued regions and regions of likely to be caustic 
for this potential in Minkowski space time as well as in FRW expanding universe which broadly agrees with the analysis of 
Ref. 1I15I1 . However in the case of the expanding universe, the caustic formations is more definite than the case of 1 H- 1 
dimensional Minkowski space. It is clearly seen from figure |4] that expansion works against caustic formation: it dilutes the 
effect of caustics but can not render the situation caustic free in the present case which is generally true for a tachyon potential 
that decays faster than at infinity. 

The field exhibits oscillatory behavior for exponentially increasing rolling massive scalar field potential, V = Vbe^^^ . 
In the case of inhomogeneous field, the magnitude of time oscillations of the field is very small and hence the field can be 
considered as almost steady for a sufficiently tuned values of initial parameters. The field with exponentially increasing rolling 
massive potential is free from caustic formation in both cases of 1 H- 1 dimensional Minkowski space and expanding universe, 
see Fig. 9. It is interesting to note that this potential can give rise to late time acceleration provided that we fine tune the energy 
scale Vq in the potential appropriately. 

For inverse power-law potentials, V{(j)) ~ Vo/4>^\ < Ji < 2, dark energy is a late time attractor of dynamics and in this 
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FIG. 12: The particle trajectories in 1 + 1 Minkowski space obtain from equation l ll6t with potential i35t and 4>i{q) = exp{—q^) corresponding 
to the solutions of the equation J37t with (/)o = 1 as shown in the middle panel of the figure [Tol Caustics and multi-valued regions are clearly 
formed in the space. 



case, we did expect the system to be free from caustics. Our analysis shows that the time variation of the homogeneous and 
inhomogeneous fields are almost similar in this case, because the space variation of the field is insignificant in comparison to 
time variation. Fig. 12 shows that caustics are formed with the multi-valued regions beyond them in the field configuration in 
case of inverse power-law potentials in Minkowski space time. 

Heuristically speaking, expansion generally works against caustic formation and its effects becomes crucial for inverse power- 
law potentials with < n < 2. If the potential vanishes faster than l/c/)^ at infinity, the dust like solution is late time attractor; 
the exponentially decreasing potential belongs to this category. In this case, caustics form in Minkovski space very distinctly. 
The cosmic expansion does dilute the effect of caustic formation but can not irradicate them. It is really interesting that for the 
inverse power law potentials under consideration, the effect of expansion can compete with the tendency of caustic formation. 
Indeed, dark energy as a late time attractor of the dynamics in this case, gives rise to cosmic repulsion allowing to avoid caustics 
and multi-valued regions in the field profile. 

From the behavior of the scalar field with exponentially increasing and inverse power-law potentials, we infer that field first 
rolls fast and mimics dark matter and subsequently gives rise to dark energy, see figure [T3] Our simulation shows that the 
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FIG. 13: Numerical solutions of equation il9\ for the three different field potentials (left). The dotted line in the figure indicates the infinite 
value. The time variations of the equation state parameter w for all three potentials (right). These plots are obtained by taking constants 
k = 8nG, Vo, M, n as unity and the initial field value 00 = 0.1. 




FIG. 14: The particle trajectories in FRW expanding universe obtain from equation J23b assuming all time dark energy effect, i. e. with the 
scale factor a ~ t2/3(i+«,) ^ _ Q g ^^^^^^ 

using the numerical data of solutions of equations (18) and (19) (right) for the inverse 
power-law field potential i35i and the profile (pi (q) = exp{—q^). Other assumptions for the right panel of this this figure are same as the right 
panel of the figure [9] We put both left and right panels in this figure to focus on the late time evolution of the field as well as to show that at 
late times the inverse power-law potentials have exactly same effect on the scalar field as the dark energy has on it. 



particle trajectories computed in the expanding universe using the exact evolution equations (18) and (19) broadly agree with 
result obtained assuming dark energy dominance at all times. 

Evolution of field configurations at different times is similar for the cases of exponentially decreasing rolling massive and 
inverse power-law potentials apart from the initial and transition periods in the Minkowski space and the expanding universe. 
It should be noted that in the case of the expanding universe, the pattern of the evolution of the field remains almost same as 
its initial profile for the scalar field dominated universe. Since the field is almost steady for the exponentially increasing rolling 
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massive scalar potential, the field configuration is same as its initial profile. Thus taking into account the results of Ref. ifisll . 
we conclude that caustics generally form for DBl systems in Minkowski space time with an exception of massive rolling scalar 
On the other hand in FRW expanding universe, caustics form in DBl systems with potentials decaying faster than 1/0^ at 
infinity, in particular the exponentially decreasing rolling massive scalar potential. As for the inverse power-law potentials under 
consideration, they are free from caustics and may be suitable to explain the late time cosmic acceleration. 
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